Nonsingular Increasing Gravitational Potential for the Brane in 6D 
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Abstract 

We present a new (l+3)-brane solution to Einstein equations in (l+5)-space. As 
£N| , distinct from previous models this solution is free of singularities in the full 6-dimensional 

space-time. The gravitational potential transverse to the brane is an increasing (but 
not exponentially) function and asymptotically approaches a finite value. The solution 
localizes the zero modes of all kinds of matter fields and Newtonian gravity on the brane. 
An essential feature of the model is that different kind of matter fields have different 
localization radii. 



PACS: ll.10.Kk, 04.50.+h, 98.80.Cq 



The scenario where our world is associated with a brane embedded in a higher dimensional space-time 
with non-factorizable geometry has attracted a lot of interest since the appearance of papers |3 El ■ 
| In this model gravitons, which are allowed to propagate in the bulk, are confined on the brane because 

of a warped geometry. However, there are difficulties with the choice of a natural trapping mechanism 
for some matter fields. For example, in the existing (l+4)-dimensional models spin and spin 2 fields 
are localized on the brane with an exponentially decreasing gravitational warp factor, spin 1/2 field are 
localized on the brane with an increasing factor 0], and spin 1 fields are not localized at all 5 . For the 
case of (l+5)-dimensions it was found that spin 0, spin 1 and spin 2 fields are localized on the brane with 
a decreasing warp factor and spin 1/2 fields again are localized with an increasing factor UJ. 

The reason why there are problems with localization of fermions in warped geometries is that in the 
Lagrangians for the fermions there appears an increasing exponential, coming from the metric tensors with 
upper indices, g AB , and from the tetrads, h^. As a result the action integral over the extra coordinates 
diverges, which is the signal for the non-localization of the fermionic fields. In both (1+4)-, or (1+5)- 
space models with warped geometry one is required to introduce some non-gravitational interaction in 
order to localize all the Standard Model particles. 

For reasons of economy and to avoid charge universality obstruction [7] one would like to have a 
universal gravitational trapping mechanism for all fields. In our previous papers we found such a solution 
of the 6-dimensional Einstein equations in (2+4)- and (l+5)-spaces, which localized all kind of bulk 
fields on the brane These solutions contain non-exponential scale factors, which increase from the 

brane, and asymptotically approach a finite value at infinity. In the paper |1(J| the solution of |§] was 
generalized to the case of n dimensions. In this paper we present a new solution to the Einstein equations 
in (l+5)-space, which, similar to the models [SJ|5] , also localizes all kind of physical fields on the brane, 
but is free of singularities in the full space-time. Because of this feature of the solution by setting realistic 
boundary conditions we are able to fix the free parameters of the model. 
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The general form of action of the gravitating system in six dimensions is 



d 6 x\f ~ u g 



-6, 



M 4 

— ( 6 i? + 2A) + 6 L 



(1) 



where a/ — 6 g is the determinant, M is the fundamental scale, 6 R is the scalar curvature, A is the cosmo- 
logical constant and 6 L is the Lagrangian of matter fields. All of these quantities are six dimensional. 
The 6-dimensional Einstein equations with stress-energy tensor Tab are 

6 Rab - \ 9 ab *R=-^Z (^9ab + Tab) ■ (2) 

Capital Latin indices run over A, B, ... = 0, 1, 2, 3, 5, 6. 

As in the papers |SJ |U| for the metric of the 6-dimensional space-time we choose the ansatz 

ds 2 = <j) 2 (r)T la/3 (x v )dx a dx 13 - X(r)(dr 2 + r 2 d9 2 ) , (3) 

where the Greek indices a, /?,... = 0,1,2,3 refer to 4-dimensional coordinates. The metric of ordinary 
4-space, rj a p{x v ), has the signature (+, — , — , — ). The functions <j>(r) and A(r) depend only on the extra 
radial coordinate, r, and thus are cylindrically symmetric in the transverse polar coordinates (0 < r < oo, 
< < 2tt). 

The ansatz (0 is different from the metric investigated in (l+5)-space brane models with warped 
geometry |51 ITT1IH] 

ds 2 =4> 2 {r)T la0 {x v )dx a dx fs ~dr 2 -\{r)d6 2 . (4) 

In (pj) the independent metric function of the extra space, A(r), serves as a conformal factor for the Eu- 
clidean 2-dimensional metric of the transverse space, just as the function <j> 2 (r) does for the 4-dimensional 
part. However, in J3J the function A(r) multiples only the angular part of the metric and corresponds to 
a cone-like geometry of a string-like defect with a singularity on the brane at r — 0. Only in the trivial 
case A = r 2 is the metric (@J regular on the brane. 

The stress-energy tensor Tab is assumed to have the form 

Tuv = -9iiuF(r), = -gijK(r), T ifl =0. (5) 

Using the ansatz (|3J|, the energy- momentum conservation equation 

V A T AB = -^=d A {^gT AB ) + Tq D T cd = , (6) 
V 9 

gives a relationship between the two source functions F(r) and K(r) from JSJ 

K' + 4^- (K - F) = . (7) 

We want to point out a problem associated with the source functions (0). In general the Einstein 
equations have an infinite number of solutions generated by different matter energy-momentum tensors, 
most of which have no clear physical meaning. There is a great freedom in the choice of F(r) and 
K(r); there are no other restrictions, except Q, on their form. It is not easy to construct realistic 
source functions from fundamental matter fields so that the brane is a stable, localized object. We shall 
determine the sources F(r) and K(r) from some general physical assumptions that they are smooth 
functions of the radial coordinate r, describe a continuous matter distribution for all r, and that they 
decrease outside the brane r > e, where e is the brane width. 

For a string- like defect, which corresponds to the metric ansatz l|4|l. a set of n scalar functions with a 
Higgs potential (which breaks the global SO(n) symmetry to SO(n — 1)) are used as a source to make the 
brane stable. In this case topological arguments guarantee the stability of the brane because the (n — 1) 
homotopy group is the integers: II„_ 1 (50(n — 1)) = Z (see, for example [§])■ 
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In our case we require the condition of classical stability |13| , that the total momentum of the brane- 
matter configuration in the direction of the extra dimensions is zero 



P l = J t?dS A = , (8) 



where is the total energy-momentum pseudo-tensor of gravitation plus matter fields on the brane. For 

B 
A 



t~A one can choose, for example, the so-called Lorenz energy-momentum complex 



^{g BD g CE (d D g AE -d E g AD )} . (9) 



To obey the stability condition (JHJl, brane solutions must satisfy the conditions 

t? = i = o . (io) 

Using JHJ from l|10(l we obtain that the metric tensor of the transverse space is a function of r only and 
that giA = 0. These conditions are satisfied by both ansatze © and (@J. One can have a condition 
similar to H1(J|) separately on the gravitational energy-momentum pseudo-tensor, which gives ^3] 

W% = , (11) 

where are components of the 6-dimensional Christoffel symbols. From follows the standard form 
of the 4-dimensional part of the brane metric tensor, g a p = (j) 2 (r)ri a 0(x v ). For the case of more then six 
dimensions 9 dependent coefficients will appear in the angular part of the metrics © and (@J) and the 
condition of classical stability (JSJ) fails. This may be one of the reasons why brane solutions in more than 
six dimensions are singular [T%| . 

The condition of classical stability JSJ) is satisfied for both metrics Q and |0J. It is expected that, as 
for the string- like defect Q), one can find for particle physics models for the source (JHJ corresponding 
to stable brane configurations. 

To solve equations (O we require that the 4-dimensional Einstein equations have the ordinary form 
without a cosmological term 

- \ih»R = o • (12) 

The Ricci tensor in four dimensions R a p is constructed from the 4-dimensional metric tensor rj a fj{x v ) in 
the standard way. Then with the ansatze © and 10 the Einstein field equations (0) become 

„<!>" , 1A" 1(A') 2 1 A' A 

f + ^ + ^ = ^>-M. 

where the prime = d/dr. These equations are for the aa, rr, and 99 components respectively. 
Subtracting the rr from the 99 equation and multiplying by </>/</>' we arrive at 

This equation has the solution 

A(r) = i± , (15) 
r 

where p is an integration constant with units of length. 

System (|13|) . after the insertion of l|15|) . reduces to only one independent equation. Taking either the 
rr, or 99 component of these equations and multiplying it by r0 4 gives 

+ fa' + 3r0 2 (0') 2 = P ^r[ K ir) A] . (16) 
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In our previous paper jS] the source functions F(r) and K(r) outside the core r > e were taken to 
have the form 

F(r > e) = K(r > e) = L , (17) 
where / is some constant. Then from the Einstein equations the following solution was found 

r b - c b 

<^ = V+^' (18) 

where 

[W , aAp 2 b b a-l 

a = \ — - , o = 7 , c — e (19) 

V 3A 5M 4 ' a + 1 y ' 

are integration constants. The solution l|18|) is an increasing function from the brane to some finite value 
at infinity 

#oo) = a = ^ > 1 . (20) 

The factor l/<fi 2 (r) has 5-\ike behavior outside the core and the source functions (|17fl decrease as required. 

In P| it was shown that the solution l|18H provides a universal, gravitational trapping for all kinds of 
matter fields. However, in this model we did not specify source functions on the brane < r < e and 
there were a large number of free parameters. Now we want to choose the source functions F(r) and 
K(r) everywhere, so that the solution (j> will localize all kind of physical fields on the brane and be a 
regular function in the full 6-dimensional space-time. 

We require for <f> the following boundary conditions near the origin r = 

<j){r -> 0) w 1 + dr 2 , <t>'(r -> 0) w 2dr , (21) 

where d is some constant. At infinity we want <j)(r) to behave as 

^(r-^oo)^a, 4>'(r oo) , (22) 

where a > 1 is some constant. Since the function <fi' is proportional to the metric of the extra 2-space, 
the boundary conditions (|22|l imply that at infinity A — ► and the effective geometry is 4-dimensional. 

The source functions F(r) and K(r), which satisfy restriction JJJ and give a desirable solution were 
found recently in the paper JS] 



where f\,f% are constants. Note that these source functions do not have a vanishing value at r — > oo, 
due to the asymptotic behavior of </> given in l|22|l . 

Substituting 123|) into i|16|) , taking its first integral and setting the integration constant to zero yields 

[HI EH] 

r ^ = 777774 TT + TT^ ~ 4> ) ■ (24) 



lOAf 4 V 3A 4A 
By introducing the parameters A and a such that 

10M 4 'J 1 5 ' 5 V ; ' V ; 

equation l|24|l becomes 

r<j>' = A[-a+ (a+ l)4>- 4> 2 ] ■ (26) 
Equation l|26() is easy to integrate ^5] 



c b + ar b 



(27) 
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where b = A(a — 1) and c are integration constants. From the boundary conditions (|21|l it follows 



b = A(a - 1) = 2 . (28) 

The width of the brane e corresponds to the inflection point of the function </> (where the second 
derivative of <f> become zero and then changes sign). From the condition 4>"(r = e) = we can fix the 
integration constant c in (|27[) 

c 2 = 3e 2 . (29) 

Finally the solution <j) corresponding to a non-singular transverse gravitational potential for the brane 
has the form 

3e 2 + ar 2 

^ = ^T^- (30) 

From the condition that we have a 6-dimensional Minkowski metric on the brane, \(r = 0) = 1, 
(any other value corresponds only to a re-scaling of the extra coordinates) we can fix also the integration 
constant in (|T5)l 

3f 2 

p 2 = *nrn ■ (31) 

2(a — 1) 

Then using (|28|l the brane width can be expressed in terms of the bulk cosmological constant and funda- 
mental scale 

Now the metric tensor of the transverse space (|15|l is not dependent on a and has the form 

9e 4 

A = ( 3e 2 + r 2)2 ' (33) 

Using solutions H27J) . I|33|l and the relationship (|15(l to integrate the gravitational part of the action 
integral QJ over the extra coordinates we find 



All f A//-4 ft* rOO p 

S g = —J dx 6 ^f^g 6 R=— J dd J dr r(j> 2 \ J dx^yf^R 

/ a r fufi r 

d(f> cf> 2 J dx 4 y^fjR = — e 2 7r(a 2 + a + 1) J dx 4 y^jR , (34) 

where R and rj are respectively the scalar curvature and determinant, in four dimensions. 

The formula for the effective Planck scale in our model, which is two times the numerical factor in 
front of the last integral in l|34|l 

m 2 Pl 

is similar to those from the "large" extra dimensions model . The differences are, the presence of the 
value of gravitational potential at extra infinity, a, in (|35|l . and that the radius of the extra dimensions 
is replaced by the brane width e, which, as seen from l|32() , is expressed by the ratio of the fundamental 
scale M and the cosmological constant A. 

The normalization condition for a physical field, that its action integral over the extra coordinates r, 6 
converges, is also the condition for its localization. As was shown in |2j Newtonian gravity is localized 
on the brane, since the action integral for gravity, 134|) . is convergent over the extra space. However, the 
wave-functions of a localized matter field can be spread out from the brane more widely then the brane 
width e. In order not to have contradictions with experimental facts, such as charge conservation [7], the 
parameters of the model must be chosen in a proper way. 

When wave-functions of matter fields in six dimensions are peaked near the brane in the transverse 
dimensions there wave-functions on the brane can be factorized as 



M 4 ire 2 {a 2 + a + 1) , (35) 



E(x A ) = ^-L , (36) 
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where the parameter re is the value of the constant zero mode with the dimension of length. These 
parameters can be found from the normalization condition for zero modes 



2tt 



1 



d9 / dr^g- = yf=j , (37) 



which also guarantees the validity of the equivalence principle for different kinds of particles. 

Let us consider the situation with the localization of particular matter fields. If we assume that the 
zero mode of a spin-0 field, in six dimensions is independent of the extra coordinates its action can be 
brought to the form [5] 

S* = J tfxyf^g & L^>{x A ) = J dr r(f> 2 \ J d^x^L^x") = 

= ^! f 4> 2 d4 f d 4 xV=vL*(x u ) = ± a + 1} f d^x^L^x") , (38) 

K $ J\ J K $ J 

where L${x v ) is the ordinary 4-dimensional Lagrangian of the spin-0 field and is value of the constant 
zero mode. The integral over r, 6 in (|38|l is finite and the spin-0 field is localized on the brane. 

The action for a vector field in the case of constant extra components (Ai = const) also reduces to 
the 4-dimensional Yang-Mills action multiplied an integral over the extra coordinates 9 

S A = J dPxy/^g 6 L A (x B ) = J o dr r\ J d A x^L A {x v ) 

J"d<f> j d A x^)L A {x v ) = ^ J d 4 x^L A {x v ) , (39) 



2irp 2 



where k a is the value of the zero mode of the vector field. The extra integral in Ij39(l is also finite and 
the gauge field is localized on the brane. 

The factorization of the zero mode of a 6-dimensional spinor field in the ansatz © is different from 
the definition (|36|) . having instead the form [5] 

HxA) = - (4o) 

(r<p ') 

where is the value of the constant zero mode. Integrating the 6-dimensional action of fermions over 
the extra coordinates, using the explicit form (|30|l . yields 



S* = j d e xy/^ 6 L^(x A ) = ^ jf dr J d 4 x^jL*(x») 

d^x^r\U ! {x v ) , (41) 



3tt 2 6 2 



K%y/2a(a - 1) 

where is the 4-dimensional Dirac Lagrangian. The extra 1/0 dependence in the second integral of 
(|41|) comes from the tetrad functions with upper index in the definition of the Dirac gamma matrices for 
the ansatz @. The integral in (|41|l over r and 6 is finite and Dirac fermions are localized on the brane. 

Equating the coefficients of action integrals ||3SJ|, l|39|) and l|41f) to 1, so as to satisfy the normalization 
condition i|37|) , and to guarantee the equivalence principle for gravity, we find the values of the zero modes 
for spin 0, spin 1 and spin 1/2 fields 

3tt 2 € 2 

K% = ire 2 (a 2 +a + 1) , n\ = ine 2 , n% = , (42) 

y/2a(a — 1) 

which are used to parameterize the 4-dimensional fields in the Lagrangians. 

Within our model we now want to find the the positions where the zero modes are localized as well 
as the localization radii of the different fields. 
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From 138|) the effective zero mode wave-function of the scalar field in flat space can be defined as 



$o(r) 



'2irr(j) 2 \ 



27r3e 2 



: (3 £ 2 



(3e 2 



(43) 



where k$ has the value l|42|) . Function l|43|l is zero at infinity (r — > oo) and on the brane (r = 0) and has 
a maximal value at some localization distance, d$, between the brane (r = 0) and infinity (r = oo). This 
localization distance, d$, and the localization radius r$, can be found by equating the first and second 
derivatives (to find the maximum and inflection point of the function) of l|43|l to zero respectively. For 
the localization distance this yields 



, 5a - 7 + V49 - 58a + 25a 2 
2a 



(44) 



From this formula we see that since a > 1 the maximum of the wave-function of scalar fields 1431) is 
located outside the brane d$ > e. 

Setting the second derivative of l|43|) to zero gives 



5ar 6 + (63 - 66a)eV - (102 - 45a) eV 2 - 9e 6 







(45) 



This is an effectively cubic equation for r 2 , which has one real and two complex solutions. The radius, 
r$ , of the zero mode scalar wave-function is given by the real solution to (|45|) . The expressions of the 
solutions of l|45|). which can be obtained using a symbolic mathematics program such as Maihematica, 
are extremely long and we do not write them here explicitly. 

From l)39|l the effective wave-function of the vector field zero mode takes the form 



Mr) 



'27rrA 



27r3e 2 



ha 3e 2 



(46) 



This function is also zero on the brane and at the infinity, and has a maximal value at some distance, 
cLa, in between. Again setting the first and second derivatives of l|4^|) to zero we find 



I'A 



e 2 (9 + 4V6) 



1.9e 



(47) 



So the peaks of the vector wave-functions are located exactly at the edge of the brane, r = e and the 
radius of localization is approximately 2e. 

For the effective fermionic zero modes from (|41ll we have 



M r ) 



'2tt P 2 



rcj) 2 



1/2 



1/4 



1 



V3i 



(48) 



This function is zero at infinity, but unlike the wave-functions of the scalar and vector zero modes, the 
peak of the fermion wave-function coincides with brane location, r = 0. From the inflection point of the 
function l]48[) . which is found by equating the second derivative of (|48|l to zero, we obtain the localization 
radius for fermions 

r* = ey/I. (49) 

From this formula we see that since a > 1 the peaks of fermionic wave-functions on the brane are very 
sharp. 

After introducing interactions between the various fields one might be able to use the various overlap- 
pings of the wave functions of the spin 0, spin 1 and spin 1/2 fields in speculations to explain the different 
types of mass hierarchies in particle physics. For example, in the split fermion model Hi localization of 
different species of fermions (not fields with different spins as in our case) at different points of a single 
thick brane was used to solve the hierarchy problem. 
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To summarize, in this paper it is shown that for a realistic form of the brane stress-energy, there exists 
a static, non-singular solution of the 6-dimensional Einstein equations, which provides a gravitational 
trapping of 4-dimensional gravity and matter fields on the brane. An essential feature of the model is 
that different kinds of matter fields have different localization distances from the brane. This property is 
in principle experimentally testable. 
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